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Abstract 

It is shown using a method based on a modified version of the 
mean field theory of Miklavic Marcelja that it should be possible for 
osmotic pressure due to the counterions associated with the two poly- 
electrolyte polymer brush coated surfaces to support a reasonable load 
(i.e., about 10 6 Pa) with the brushes held sufficiently far apart to pre- 
vent entanglement of polymers belonging to the two brushes, thus 
avoiding what is believed to be the dominant mechanisms for static 
and dry friction. 

1 Introduction 

Polymer brush coatings on solid surfaces provide very effective lubrication, 
in the sense that they are able to support significant load (pushing the sur- 
faces together), but have exceedingly low friction coefficients [lj. Human and 
animal joints are known to exhibit very low friction and wear. The outer 
surface of the cartilage coating these joints have polymeric molecules pro- 
truding from them |2l El [H El El El [H El HQl HH H21 H31 [H] . This suggests the 
strong possibility that their very effective lubrication is a result of polymer 
brush lubrication. It has been shown that at small loads, polymer brush 



coated surfaces can slide relative to each other with the bulk of the brushes 
not in contact [T5]. If the surfaces are far enough apart under such loads, 
the load will be supported almost entirely by osmotic pressure due to a di- 
lute concentration of polymers that protrude into the thin interface region 
separating the brushes. As a consequence, there will be little entanglement 
of the polymers belonging to the two brushes (i.e., penetration of polymers 
belonging to one brush into the second brush). It was argued in Ref. [TH] 
that such entanglement of polymer brushes leads to static and kinetic friction 
that saturates at a nonzero value in the limit of zero sliding velocity (i.e., dry 
friction) of polymer brush coated surfaces. It is reasonable to expect that 
there should be little wear as well when there is little entanglement. 

The discussion in Ref. [T5] deals with static and slow speed kinetic fric- 
tion, and not the purely viscous kinetic friction that occurs at high sliding 
speeds. Since the polymers hyaloronan or lubricin, which coat the cartilage in 
human and animal joints [16], are charged, it is necessary to consider polyelec- 
trolyte brushes, whose equilibrium properties have been studied using mean 
field theory by Zhulina, et. al.[T7j, Misra[18j and Miklavic[19j. and Pincus 
and Tamashiro, et. aLj^O]- Raviv, et. al.|21j. have found that polyelectrolyte 
brushes exhibit remarkably low friction coefficients (10~ 3 or less) compared to 
the friction coefficients typically found for neutral polymer brushes [Tl |22| [23] . 
It is proposed in the present article that polyelectrolyte polymer brushes are 
more effective lubricants than neutral polymer brushes discussed above be- 
cause, as will be shown, counterion osmotic pressure present in a relatively 
thin interface region separating two polymer brushes is able to support a load 
of the order of 10 6 Pa (which is comparable to the loads supported by the 
polyelectrolyte brushes studied in Ref. [2TJ), without the tops of the bulk 
part of the density profiles of the brushes being in contact, as illustrated in 
Fig. 1. The height of a neutral polymer brush is determined by a competition 
between the polymer's elasticity (of entropic origin) and the mutual repulsion 
that occurs in a good solvent between the monomers making up the poly- 
mers of the brush. For polyelectrolyte brushes, the counterion entropy [20] 
also plays a significant role in stretching the polymers in the brush. 

The present treatment of lubrication due to polyelectrolyte brushes is 
based on the theory of Miklavic and Marcelja[19], which is an extension 
of the analytic mean field theory due to Milner, et. al.|24J. As has been 
shown in several treatments of the subject [25], using numerical treatments 
of mean field theory and molecular dynamics[26j, the analytic treatment of 
mean field theory [2"1] is only accurate for relatively stiff polymer brushes, for 
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Figure 1: The geometry of two polyelectrolyte polymer brush coated surfaces 
with the load pushing the surfaces together supported by osmotic pressure 
due to counterions in the interface regions separating the tops of the brushes 
is illustrated schematically. The dots located among the polymer chains and 
in the interface region between the two brushes represent the counterions. 
As illustrated here, D denotes the spacing of the surfaces and h denotes the 
polymer brush height. 
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which the brush height is large compared to the unswollen single polymer 
radius of gyration iV 1//2 a, where a is the monomer length. Since for good 
solvents this is a useful limit to consider, a treatment based on analytic 
mean field theory is expected to give accurate results in the important case 
of relatively dense polymer brushes with relatively strong repulsions between 
their monomers, as well as for polyelectrolytes for which counterion osmotic 
pressure within the brushes provides strong forces to keep the polymers in 
the brush well stretched [20]. Furthermore, analytic mean field theory shows 
important trends as the brushes are compressed [15], which are consistent 
with experiment. Therefore, a theoretical study of lubrication by polymer 
brushes based on mean field theory is expected to be a good starting point 
in an effort to understand lubricating properties of polymer brushes. In 
section 2, it is argued that the method of estimating the mean distance that 
a polymer belonging to one of the brushes is able to extend into the interface 
region separating the two brushes, used for neutral polymer brushes, applies 
equally well to polyelectrolyte brushes, when they are sufficiently compressed, 
and this distance is estimated. If the interface region separating the tops of 
the brushes is wider than this distance, polymers from one brush will be 
unlikely to penetrate into the second. In section 3, the question of how much 
load can be supported by counterion osmotic pressure (with the brushes held 
sufficiently far apart to prevent the occurrence of static friction) is studied. 

In an effort to determine a possible mechanism for the ultra-low friction 
found in Ref . [2TJ , an approximate solution for the Poisson-Boltzmann equa- 
tion beyond the Debye-Huckel approximation will be used to determine the 
concentration of counterions in a region located midway between the two 
polyelectrolyte brushes. This result is used to show that for polyelectrolyte 
brushes, osmotic pressure due to the counterions is capable of supporting 
about 10 6 Pa of load, even with the tops of the mean field theoretic monomer 
profiles of the two brushes about 14A° apart. In fact, it will be shown that 
this result is valid even for salt concentrations smaller than but still com- 
parable to 0.1 M (about 10 26 m~ 3 ), the salt concentration in living matter. 
Increasing the salt concentration further, however, will put the system in a 
regime in which Debye-Huckel (D-H) approximation is accurate and predicts 
that the ion concentration mid-way between the brushes is negligibly small 
once the D-H screening length is much smaller than the width of the interface 
region. 
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2 Estimate of the Brush spacing for which 
Static Friction due to Brush Entanglement 
will Not Occur 

Part of Ref. [TS] relevant to the present discussion is briefly summarized 
in this paragraph. Milner, et. al.[24j, proposed a simple way to solve the 
mean field theory analytically for neutral polymers attached at one end to a 
surface. In this treatment, the location of the n th monomer of the i th polymer 
belonging to a polymer brush of minimum free energy, i"i(n) satisfies the 
differential equation 

where V(ri(n)) = w0(ri(n)), and where for neutral brushes, 0(r) is the 
monomer number density and w is the strength of the monomer- monomer 
repulsion parameter. This can be thought of as an "equation of motion" for 
the monomers, in which the index n labeling the monomers formally plays 
the role of time. It is formally analogous to Newton's second law for motion 
of a particle in a potential equal to V(r(n)). Existence of a solution of the 
mean field equations of motion requires that V(r(n)) be a parabolic function 
of z [24] . In Ref. [15] it was shown that when two polymer brushes are 
in contact with each other or nearly in contact, polymers belonging to one 
brush can penetrate into the second brush. As a consequence, there is a force 
of static friction equal to the force needed to pull these intertangled poly- 
mers out. The force of static friction per intertangled polymer was argued 
to be equal to fceT/^, where £ is the mesh size of one of the two polymer 
brushes. The magnitude of the friction was estimated in Ref. [TH] to be of 
the order of 10 3 Pa or more. Even when the applied force is below this value, 
the surfaces will not remain truly stationary, but rather, will creep relative 
to each other. The reason for this behavior is that the intertangled polymers 
will diffuse out in a diffusion or reptation time r. It was shown that at high 
compression r can be sufficiently long to consider this to be a true force of 
static friction. If a force greater than this static friction force is applied, and 
hence, the surfaces slide with a speed much greater than this creep speed, 
since the polymers will no longer have enough time to re-entangle in the 
second brush, the friction will no longer be determined by this mechanism. 
Instead, there will be viscous friction resulting from the fact that the sol- 
vent gets sheared as the surfaces slide. This viscous friction force per unit 
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surface area in the slow sliding speed limit is given by rj(v/£ p ), where r\ is 
the viscosity of the solvent, v is the sliding velocity[27] and £ p is the hy- 
drodynamic penetration length into the brushes. [28]. Let us consider sliding 
speeds that are sufficiently slow so that the zero sliding speed configuration 
of the polymers is not significantly disturbed. For uncompressed polymer 
brushes, which have a parabolic density profile [21], £ p is comparable to the 
equilibrium brush height [25], but for brushes which are compressed because 
they are supporting a load, the density profile will get flattened out, and as 
a consequence, £ p can be comparable to the polymer spacing at the surface 
to which they are attached, s. Then, for example, assuming the solvent to be 
water, for v=l mm/s, the largest speed used in the experiment, the force of 
friction per unit area for the uncompressed brush case will be approximately 
equal to about a Pa, for brushes of equilibrium height of the order of 500^4°. 
For compressed brushes with an anchor spacing s of 84v4°, it will be about 
5 times this value. This estimate is an upper bound on the shear stress, 
since it assumes no-slip boundary conditions for the fluid at the polymers 
and the surface. For plate spacing just above the maximum separation for 
which the friction observed is above the experimental accuracy (reported in 
Ref. [23]), the load is reported to be F/R=0.01N/m (where R is the radius of 
the cylinders in the surface force apparatus, which is about 0.01 m), which, 
using the standard Hertz formula[29j, gives a contact area for the surface 
force apparatus of about 10~ 10 m -2 . The maximum observable shear force in 
this experiment is 0.25/iiV, which when divided by the above contact area 
gives a shear stress from the experimental data of about 10 4 Pa, clearly well 
above the shear stress of about a Pa found above for viscous friction. Hence 
the viscous friction discussed above is clearly well below the experimentally 
observed kinetic friction. In contrast, Eq. (16) in Ref. [U] shows that it 
is easy to get a value for the static friction or slow speed kinetic friction, 
due to the polymer blob entanglement mechanism discussed earlier, compa- 
rable or greater than the experimental value if the brushes are sufficiently 
compressed. 

An important advantage of polyelectrolyte over neutral polymer brushes 
as lubricants is that the osmotic pressure due to counterions might, under 
the right conditions be able to support the load, allowing the brush coated 
surfaces to slide without the bulk of the brushes being in actual contact, 
and hence, with negligible friction due to entanglement of polymers from 
one brush in the second. In order to calculate the separation of the two 
brushes above which friction due to this entanglement no longer occurs and 
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to calculate the contribution to the repulsion of two polyelectrolyte brushes 
due to osmotic pressure due to the counterions at this separation, we will 
use an analytic mean field theory treatment of polyelectrolyte brushes based 
on the treatment due to Miklavic and Marcelja[19]. In their treatment they 
use the mean field treatment of Ref. [23] with V(r) = w<p(r) + ei/)(r), where 
•0(r) is the electrostatic potential, due to the charged polymers making up 
the brush, screened by the counterions divided by the electronic charge e, 
which satisfies Poisson's equation 

^ = 4vrp(^)/e-(4vr//e)0(z), (2) 

where p(z) is the ionic charge density, z is the distance from the lower surface, 
4>(z) is the monomer density profile of the brush (or brushes, if there are two 
of them in contact or nearly in contact) and f is the charge in units of the 
electron charge e per monomer. In contrast to Ref. [19] , in the present treat- 
ment of this problem, the Debye-Huckel approximation will not be used when 
estimates of the osmotic pressure due to counterions are made. The charge 
per monomer f has a maximum value because of Manning condensation [30J. 
If the charge spacing on the polymer is less than the Bjerrum length, there 
will be Manning condensation |30j, which means that some of the counteri- 
ons will condense onto the polymer, until its charge density is reduced to 
the point that the charge spacing becomes equal to a Bjerrum length for 
monovalent counterions. Then, since the Bjerrum length, £b = e 2 /{ek,BT) 
is 7.1 Angstroms for a solvent with a dielectric constant comparable to that 
of water, the largest possible value of f is the ratio of a monomer spacing 
to a Bjerrum length, or 0.35, if we assume the value for the monomer spac- 
ing from Ref. [21] of 2.bA° . Since we want the polymer brushes to behave 
as polymer brushes in a good solvent, we want the charge per monomer to 
be sufficiently small so that there are not too many counterions condensed 
on the polymers, in order to prevent possible collapse of the brushes due to 
interaction of dipole moments resulting from counterions condensed on the 
charged monomers|32]. Milner, et. al. [23] and Miklavic and Marcelja|19] 
show that in mean field theory the height of the n th monomer on a polymer 
in a brush, z n satisfies the differential equation 

d 2 z n /dn 2 = (a 2 /k B T)(dV(z n )/dz n ), (3a) 

where for polyelectrolyte brushes V(z n ) is the potential acting on the n th 
monomer, due in the present case to both hard core intermonomer forces 
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and electrostatic forces. In order to have a self-consistent solution to Eq. 
(3a), V(z) must have the form A — Bz 2 , where A and B are constants |24j. 

Let us now estimate effects of fluctuations from mean field theory by a 
method similar to the methods of Ref. [21]. If ji monomers belonging to 
a polymer are either pulled out or thermally fluctuate out of a brush, the 
solution for function of n becomes for z n < h 

z n = (h/cos(ujji))cos(ujTi), (36) 

using Eq. (5b) of Ref. [15] and requiring that Zj x = h, where j\ is the value 
of n such that for n < ji, z n > fe[l5j. For z n > h the "force" in the "equation 
of motion" is no longer determined by —(^4 — Bz 2 ), but by the electrostatic 
potential ip(z) outside of the brush, which is a smooth function of z. Then 
for n < ji, z n is given by z n = f(n), where f(n) represents the solution to the 
"equation of motion" with the potential given by the value of —ip(z) outside 
the brush instead of — (A — Bz 2 ). If j% << N (a condition for the validity 
of mean field theory), we may expand f(n) to lowest order in n for n < j\. 
Hence, 

h = z n = f( 3l ) « f(0) + /'(0)ja = z + f'(0)j u 

or 

z = h-f'(0)j 1 . (4) 

The assumption made here that Z n IS cl slowly varying function of n in the 
region outside the brushes is valid since the distance over which the electro- 
static potential varies is of the order of the thickness of the interface region 
separating the brushes and the distance that the polymers belonging to a 
brush stick out into this region is smaller than this distance, because we have 
chosen to consider the situation in which the spacing between the brushes is 
greater than the distance that the polymers belonging to a brush protrude 
into this region. We must require that dz n /dn be continuous at n = j\. It 
follows from Eq. (3) that 

dz n 

-^|n=o ~ /'(0) = -hutan(uii) w -hu 2 ji. (5) 

This implies that zq m h + hu 2 j 2 , and hence, dzo = 2huj 2 jidj\. As in Ref. 
|15j . we use the fact that (ksT / a 2 )(dz n / dn)\ n= Q is the tension that must be 
applied to the free end of the polymer in order to pull it out of the brush 
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to calculate the work needed to pull ji monomers of this polymer out of the 
brush, 

2 f Z ° 

AF = ((k B T/a ) I (dz n /dn)\ n=0 dz = 
Jh 

(2/3)k B T(hu 2 /a) 2 jl (6) 

From Eq. (5) we obtain the probability that ji monomers stick out of the 
polymer into the interface region, 

Poce^(-(2/3)(^ 2 /a) 2 Jl 3 ) (7) 

or using the fact that huo 2 jf = z$ — h [which is given under Eq. (5)], we find 
that 

P oc expl-C-^) 3 / 2 }, (8) 

where ^/h = (3/n) 2 ^ 3 (N 1 ^ 2 a/h) 4: ^ 3 . We expect that we will get extremely low 
friction, only if the polymers that fluctuate out of one brush do not extend so 
far that they get entangled in the second polyelectrolyte brush, because that 
will result in static friction and relatively large nonviscous kinetic friction[15j. 
For D — 2h = z^ — h » £, there will be no static friction due to entanglement 
of polymers belonging to one brush in the second. Using the parameters 
N=115, a = 2.5A° and h = 170v4° from Ref. [21J in the formula for £ under 
Eq. (8), we find that we must have zq — h > £ « 14. 4v4°. In fact, numerical 
work on neutral polymer brushes [25| [26] shows that the monomer density 
is generally quite small in the tails on the monomer density distribution of 
length £, implied by Eq. (8), and there is every reason to assume that the 
same will be true for charged polymers as well, making it likely that it is 
not necessary for Zo — h to be much larger than £ in order for the monomer 
density in the tails to be in the dilute regime [15]. in which friction due to 
blob entanglement does not occur. 

The model discussed in this section assumes that the polyelectrolyte 
brushes have a height which is a relatively small fraction of the fully stretched 
length of the polymers, as is the case for neutral polymer brushes. This would 
be the case only for polyelectrolyte brushes for which the charge on the poly- 
mers is a very small fraction of the maximum charge allowed by Manning 
theory [50]. or for polymers with a higher charge but carrying a high load, 
so that the brush heights are compressed to a small fraction of the heights 
expected for uncompressed brushes on the basis of Pincus' arguments[20j. 
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For highly charged polyelectrolyte brushes under light loads, for which the 
polymers are extended to close to their fully extended length, the high de- 
gree of entanglement needed to give the static and relatively high kinetic 
friction discussed in Ref. [15] for neutral brushes is unlikely. In fact, it was 
argued earlier that this is likely to be the reason that experimental studies 
of polymer brushes [23J show that the friction is negligible until the brushes 
are under relatively high loads. 

3 Load Carrying Ability of Polyelectrolyte Brushes 
in the Low Salt Concentration Limit 

Since in the absence of excess salt, the conterion contribution to the osmotic 
pressure falls off quite slowly with plate separation [20J, the possibility will be 
explored here that a reasonably large load can be supported by counterion 
osmotic pressure while the plates are sufficiently far apart to prevent entan- 
glement of the polymers belonging to the two brushes. Therefore, let us now 
consider the situation in which there is a sufficiently low concentration of ex- 
cess salt. In this situation, we must use solutions of the Poisson-Boltzmann 
equation beyond the Debye-Huckel approximation, since that approximation 
does not accurately describe the problem[20l [35]. Let us first consider the 
situation in which there is no excess salt, and only counterions are present. 
In the interface region (i.e., h < z < D — h), the electrostatic potential must 
be a solution of the Poisson-Boltzmann equation 

d 2 <4)/dz 2 = -Annoie^^, (9) 

where no is the counter-ion density midway between the plates, i.e., at z=D/2 
[35] . The solution to Eq. (9) for this geometry in the interface region between 
the two polymer brushes may be used to estimate the counterion contribution 
to the osmotic pressure supporting the load. The two polymer brushes may 
be formally replaced by two equally charged flat plates a distance D-2h apart 
if we impose the boundary condition that —dip/dz\ z= h = cr e as required by 
Gauss's law, where a e is the total charge (polymer charge plus counterion 
charge) contained within the brush (i.e., in the region < z < h). This 
solution is [35] 

iP = £- 1 log(cos 2 k (z-D/2)), (10) 
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where k 2 . = 2nn £ B , and the counter ion density is given by 

71 = n ^ B " = 2w W ° n/9V 
cos 2 k (z — D/2) 



Since 



D-h [D-h d z 

dzn(z) = n 



Ih cos 2 k (z - D/2) 

(2n /k )tank (D/2 - h) = 2a e , (12) 

tanko(D/2 — h) — k a e /no, (13) 

for ko<r e /n >> 1, we get the maximum possible value of k , which gives the 
maximum value of no, namely, ko = [tt/(D — 2h)]. From the definition of ko 
below Eq. (10), it follows that the largest possible value of 

"» = (i/2) MdW (14) 

which gives for the counterion contribution to the osmotic pressure 

Posm = k B Tn t = ^m £B{D B ^ 2h)2 - (15) 

For a value of the parameter D-2h, comparable to £, estimated in section 
2, or about 14.4A , we find that P osm = 2.7 x 10 6 N/m 2 , and it is inversely 
proportional to the square of the spacing between the tops of the brushes, 
i.e., the width of the interface region. 

In order to determine if the right hand side of Eq. (13) is much greater 
than one [which is the condition for the validity of Eqs. (14) and (15)], 
Eq. (2) will now be integrated, in order to determine whether a value of no 
(and corresponding value for the osmotic pressure given above for the region 
separating the tops of the brushes) comparable to the value given above will 
occur. The Poisson-Boltzmann equation [Eq. (2)] for the case of no excess 
salt can be written as [201 



j2 / 

^ = -47rf B [noe-*-M2)], (16) 

where = Ib^>- In order to make it possible to integrate Eq. (16), we 
will approximate the monomer density of a brush <f)(z) by the step function 
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4>(z) = (N/h)a8(h — z), where 9{x) = 1 for x > and for x < 0. This 
is a reasonable approximation because we are considering polymer brushes 
which are compressed because they are supporting a load, and under such 
circumstances, the parabolic density profile of the uncompressed brush gets 
flattened into a form that is not too different from the step function form 
given above[T5] [21]. Multiplying Eq. (16) by d<p/dz and integrating, we get 




8ir£ B n (e-^ z) - e~^ h) ) + {8-K£ B Naf '/h)[<f>(z) - <j>(h)] = 

8Tr£ B n (e~^ z) - 1) + 8-Ki B Naf /h$(z) - <j>{h)\, (17) 

for z between and h, using Eqs. (10-13) to simplify this expression. For 
z > h, (j>(z) = ln(cos 2 k (z — D/2), the solution described in Eqs. (10-13) [35]. 
Thus, Eq. (17) leads to the integral 



4>{z) 



m [(e-^-l)+r(0-0(/i)] 1 /2 



2k (z - h) 



z — h , . 

for z < h, where r = K^/hy, where K$ = 2ii£ B fN/(hs 2 ), the square of the 
inverse screening length within a brush, where N is the number of monomers 
in a single polymer and s is the mean spacing of the polymers of one brush 
along the surface to which they are attached and a = ko(D/2—h). Evaluation 
of this integral allows us to obtain <f)(z), if we know k . The total charge 
between the plates consists of the sum of the charge on the polyelectrolyte 
brushes and the counterion charge, is zero. Applying Gauss's law, using 
the fact that symmetry demands that the electric fields at both plates are 
equal in magnitude and opposite in direction, we find that the electric field 
at the plates is zero, if no counterions condense on the surfaces at z=0 and 
z=D. Then, we must demand that at the location of the lower plate, z=0, 
d(j)(z)/dz = 0. (It is clearly also zero at z=D by symmetry.) We can find 
4>(z = 0) by setting z=0 in Eq. (17), and setting d<j)(z) /dz\ z =o = 0. (Strictly 
speaking, this condition is only precisely correct if the surfaces belong to thick 
solids, so that there can be no solution present inside these solid blocks. If 
this were not the case, d<ft(z)/dz would not be required to vanish at precisely 
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Figure 2: The functions f \ = - 1 and /2 = r(0(L>/2 - h) - 0(0) are 

plotted versus the dimensionless potential 0(0) for r=ll. All quantities are 
dimensionless. 

z=0, but rather approximately at a short screening distance below z=0. The 
results are not expected to be modified qualitatively from what we will find 
in this section if we were to take this into account, however. In any case, this 
is not an important case for most applications, because there is normally no 
solution on the outer side of the two surfaces because there is normally solid 
material located there.) Using the value of 0(0) found in this way, we may 
then determine ko by integrating Eq. (18) with z=0. From the definitions of 
k and K under Eqs. (10) and (18), respectively, the ratio r = (K /k ) 2 is 
equal to the ratio of the unscreened charge density, fN/ (hs 2 ), of the polymer 
brush to n . A graphical solution for 0(0) is illustrated in Fig. (2) for a 
reasonably large value of r = (K /ko) 2 . We must have K /ko > 1. This is 
evident from the fact that (as is seen in Fig. 2) for K /k < 1 we cannot 
find a solution for 0(0) from Eq. (17) by setting d(j)(z)/dz\ z= o = 0. We see 
that for any value of r > 1, there exist two solutions for 0(0), one with 
0(0) close to 0(/i) and one with |0(O)| > \4>(h)\. Since in the large r limit, 
we expect the charge residing on the polymers within the polymer brushes 
to be highly screened, we do not expect 0(0) and <j)(h) to differ by very 
much. Therefore, the solution with 0(0) comparable to 4>{h) is the physically 
correct solution and the solution with |0(O)| > \4>{h)\ is rejected as being an 
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unphysical solution of the Poisson-Boltzmann equation. 

Eq. (18) can be integrated numerically. In order to facilitate this, fol- 
lowing the discussion in appendix A of the second article in Ref. [20J, we 
first integrate Eq. (18) by parts to eliminate the singularity in the integrand, 
giving 

tank (D/2- h) 
Kl/kl - 1 - tan 2 k (D/2 - h) 

m , 7 e-*\e-* -l + rU-dfh)] 1 / 2 

d( P 7 

4>{h) (r - e-vy 

Eq. (19) may also be written as 

h atana 



k h. (19) 



D/2-h (3 2 - a 2 (tan 2 a + I) 
-a- ^\^- l + ^ - m ^ 

where a = k (D/2 — h) and f3 = K (D/2 — h). Several solutions of Eq. (20) 
are given in table I. In each case, we chose the largest value of a for which 
there is a solution to Eq. (18) with z=0 for which d(f>(z)/dz\ z= o = 0. It was 
found to be difficult to find a solution for D/2 — h < h. It is clear that it 
must be possible to have solutions to Eq. (20) with D/2 — h < h. It is clear, 
however, that in order to have a solution to Eq. (18) with D/2 — h < h, 
the left hand side of this equation must be greater than 1, unless a is small, 
and we already saw for small values of a, Eq. (20) gives values of (D/2- 
h)/h which are greater than 1. For reasonably large values of r, the only 
way that this can occur is if there is a large contribution to the integral 
from (f) close to 0(0), where the integrand has an integrable singularity. Let 
us then examine the contribution to the integral in Eq. (18) for this part 
of the range of integration. To accomplish this, we write <f) = (f)(0) + <$' and 
perform the integral over a range of <$' which is small compared to <f)(h) — 0(0) . 
Then this contribution to the integral can be performed by expanding the 
function under the radical sign in the denominator of the integrand in a 
Taylor series to second order in 0', giving /(0) = — 1 + r(0 — <f)(h)) w 
d//d0|^ = ^( O \0' + (1 /2)d 2 f / d(j) 2 \^ = ^(f)' 2 . Using this, the contribution to the 
integral in Eq. (18) for close to 0(0) can be written as 



~ 1/2 Jo [0^ + V]V^ r " 1/2[fa(2( ^ )) -^ /2)] ' (21) 
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where T = (l/2)d 2 f/d<P% =m , TS = df/d<P\^ =m and S « X « <j>(h) - 
0(0). From the definition of 5, we see that 5 = e^df /d(p\^ = ^ and since 
/(</>) = fl — f2, defined in Fig. 2, we see that S becomes smaller and smaller 
as the curves fl and f2 in Fig. 2 become more and more nearly tangent 
at = 0(0). From Eq. (21), we see that we can make (D/2 — h)/h as 
small as we wish by making 5 smaller [and hence making the curves closer 
to being tangent at their point of intersection at 0(0)]. For example, for 
the value of (3 = 1.17, which we argued earlier to be appropriate for the 
polyelectrolyte brushes studied in Ref. [21], we find that the curves will be 
nearly tangent for a slightly larger than 0.65. This value of a corresponds 
to n = [O.Q5/(2n)]/[i B (D/2 - h) 2 } = 0.7 x 10 26 m~ 3 . Thus, we conclude that 
the density of counterions residing between the tops of the two brushes is 
sufficient to provide enough osmotic pressure to support a load of k B TriQ = 
0.28 x 10 6 Pa for the polyelectrolyte brushes studied in Ref. [21]. 

The degree of compression can be estimated from the discussion in Ref. 
[2D] . If the brushes are being pushed together with a compressional force or 
load per unit area P, we have by the scaling arguments of Ref. [2D] that 

fNk B T hk B T 
s 2 h Na 2 s 2 ' 

Setting P equal to the osmotic pressure in the interface region between the 
two brushes, ksTriQ, we may write the solution of Eq. (25) for h as 

h = f 2/3 Na[(l + f-\as 2 n /2) 2 ) 1/2 - as 2 n /(2f 1 / 2 )}. (23) 

If we use the following values for the parameters: no = 10 25 m -3 , a = 2.5A° 
and s = 40A°, we find that h is compressed to about 42 percent of its 
maximum value of f l ^ 2 Na 

We will see in the next section that there also exists a solution with r 
comparable to 1, if the charge density of the unscreened polymer brushes 
fN/(s 2 h) is not much larger than (n/2)/[(D/2 - h) 2 £ B }. 

4 Load Carrying Ability of Polyelectrolyte Brushes 
with Purely Electrostatic Interactions 

Zhulina and Borisov have published an analytic solution of the mean field 
theory of polyelectrolyte brushes[2S], assuming that the only interaction be- 



(22) 
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Table 1: Results of Calculations 



based on Eq. (20) 



K (D/2 - h) 


k (D/2 - h) 


(D/2-h)/h 


1.17 


0.6 


2.30 


0.70 


0.45 


1.40 


0.30 


0.35 


1.05 


1.17 


0.20 


37.7 



tween the polymers in a brush is that which results from electrostatic inter- 
action and entropic interactions resulting from the counterions. The main 
conclusion of this treatment relevant to the present work is that this model 
predicts that under such circumstances, for highly charged polyelectrolyte 
brushes practically all of the counterions lie inside the brushes. The reason 
for this is that if we assume that there are no hard core interactions between 
the monomers in the brushes (i.e., there are only electrostatic and counterion 
entropic interactions), the mean field potential is identically equal to the elec- 
trostatic potential energy of the monomers, which is equal to feip(r). Since 
the mean field potential must have the form A — Bz 2 [21], the net charge 
density in units of the electron charge e (due to both monomer and coun- 
terion charge) is given by — (e/fe 2 )d 2 ip/dz 2 = (2e/ fe 2 )B = tt 2 /(8f£ B a 2 N 2 ) : 
a constant. Since the product of this charge density and the volume of the 
brush must be equal in magnitude to the total amount of counter ion charge 
outside the brush (in units of e), for monovalent counterions, the number 
of counterions outside two brushes in contact with each other is equal to 
7i 2 hs 2 (8f£BO' 2 N 2 ). Then, clearly for highly compressed (i.e., h « f l / 2 Na 
and obviously s « f^Na) and highly charged polymers (i.e., / ~ 1), this 
model predicts that there are practically no counterions outside the brushes. 
This model predicts that for highly compressed polyelectrolyte brushes, i.e., 
brushes with heights h « f l ^ 2 Na (see Fig. 8 in Ref. [33]), the density 
profiles will approximate a step function and the tops of the two brushes will 
never be in contact (i.e. D > 2h). Hence, there will always be a thin inter- 
face region free of monomers separating the two brushes, but the thickness 
of this region will be negligibly small (i.e., less than L4°), in order to support 
a load as high as 10 5 Pa. This would imply on the basis of the arguments 
given in section 2, that highly compressed polymer brushes should exhibit 
a good deal of friction due to entanglement of blobs belonging to one brush 
in the second brush. (How the scaling arguments of Ref. [15] are modified 
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when the polymer interactions are replaced due to electrostatic, rather than 
hard core, interactions will be discussed in a future publication [33]. 

Hence, the calculations presented in the previous paragraph must give 
results consistent with Ref. [33] when (D — 2h)/D << 1. In order to make 
contact between Ref. [33J and the Poisson-Bltzmann equation calculations 
presented in the previous section, let us do such a calculation for an ex- 
ample with (D — 2h)/D << 1 with some reasonable parameters. Consider 
(K D/2) 2 = 2iTe B (2fN/(7T 1 / 2 hs 2 )D 2 1 which for £ B = 7A°, fN = W 3 A°, 
s = 10 2 A° and h « D/2 « 10 2 A°, we get K D = 15.8. The solution to 
Eq. (18) for (D - 2h)/D = 0.001 gives k D/2 = 15.79, which implies that 
n is only slightly smaller than the unscreened monomer charge density in- 
side a brush and 4>(h) — 0(0) = 0.262, which implies that <f>(z) varies by 
a relatively small fraction as z varies from to h, which is consistent with 
the calculation of Zhulina and Borisov|33j. An alternative way to view this 
is the following: In Zhulina and Borisov's solution of mean field theory [33J, 
highly charged polyelectrolyte brushes have practically no net charge inside 
the brushes, meaning that the counterion density inside the brush is nearly 
equal to the unscreened polymer charge density, as shown above. Conse- 
quently since the counterion density is continuous as we cross over from 
inside to outside the brush, k is only slightly less than K . Since (K /ko) 2 is 
equal to the ratio of n , the counterion density midway between the brushes, 
to the charge density of the unscreened polymers inside the brushes, we may 
write (-KoAo) 2 = 1 + An/n, where An is the net charge density inside the 
brushes, and n is the unscreened polymer charge density. Since An/n << 1, 
Eq. (18) reduces to 

tank (D/2 - h) « k h(An/n) (24) 

to lowest order in An/n. Since the right hand side of Eq. (18) is much less 
than 1, the tangent may be approximated by its argument, and hence we 
find from Eq. (24) that 

D/2 - h w h{An/n). (25) 

From the discussion above, we know that An = TT 2 /(8f£Bd 2 N 2 ). The un- 
screened polymer charge density is fN/(s 2 h). Using the same values for f, 
£b, a, N and h as were used above we estimate that 

D/2-h)/(D/2) ^0.4x 10" 4 . (26) 
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Thus, there exist two treatments of highly charged compressed polyelec- 
trolyte brushes, which appear to give opposite results. One predicts that 
polyelectrolyte brushes should be able to slide with their load supported 
by osmotic pressure in a region separating the tops of the brushes which is 
relatively free of monomers of sufficient thickness to prevent polymer entan- 
glement of the type that was proposed in ref. [15J to give relatively high 
friction compared to the viscous friction considered in section II and one 
which predicts that the load will be supported by osmotic pressure in a re- 
gion separating the brushes which is extremely thin. The latter would imply 
that there should always be a good deal of friction due to blob entanglement 
when h « f 1 ^ 2 Na, on the basis of Ref. [TJ]. Let us now explain this ap- 
parent contradiction. When the polymer charge is relatively high (i.e., f is 
close to 1), Zhulina and Borisov solution of mean field theory can only be 
valid if the brushes are somewhat compressed, because if the brushes are as 
stretched as far as Pincus's treatment [20J implies, the elastic free energy of 
the polymers used in Ref. [33] is not valid for highly extended polymers. 
When the brushes are compressed to heights comparable to the heights that 
they would have if the polymers were neutral, the hard core repulsion of the 
monomers will play a role comparable to that of the electrostatic interac- 
tions. Although mean field theory requires that the sum of the hard core 
interaction and electrostatic interaction potentials be equal to A — Bz 2 , with 
the constant B having the value given earlier, the electrostatic interaction 
potential alone need not have this form. In fact, when the brushes are highly 
compressed (compared to the nearly fully extended brush height bound in 
Ref. [20]) down to a height comparable to what one finds for neutral brushes, 
the hard core potential might dominate. In such a case, the electrostatic po- 
tential is certainly not constrained to have the above form. Rather, it will 
be determined completely by the solution of the Poisson-Boltzmann equa- 
tion using a monomer density profile, determined primarily by the hard core 
interaction contribution to mean field theory, in which solution like 

the one found in section III is likely to be valid. Whether or not the poly- 
electrolyte brushes are so highly compressed that hard core interaction plays 
a dominant role, the density profile of polymer brushes formed from finite 
length polymers will have significant deviations from mean field theory in 
the form of tails of length I « h at the edges of the brushes. This will 
certainly require that the mean field potential differ from the parabolic from 
required by mean field theory [21] at the edges of the brushes (i.e., near z=h) 
over a region of thickness of the order of I. Then the electric charge density 
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obtained by differentiating the electrostatic potential twice will have a charge 
density which is much larger than that found be Zhulina and Borisov[33] in 
this region, resulting in a considerably larger net charge of the polymer brush 
than that found in Ref. [33] . 



5 Effects of Excess Salt 

When there is excess salt present in the solvent, Eq. (16) gets replaced by [35] 
<h ° -47r£ B [n s (e-t - e+*) - f<f>(z)], (27) 



dz 2 



where n s represents the salt concentration and the second exponential term 
on the right hand side represents the contribution of ions with the same 
charge as the brushes to the ionic charge between the plates. For the case of 
no excess salt, described by Eq. (16), (p(z) was taken to be zero at z=D/2. 
In contrast, for the case of excess salt, described by Eq. (20), 4>(z) is usually 
taken to be zero well outside the two surfaces. For a low concentration of 
excess salt, the conditions under which Eq. (20) takes the same form as Eq. 
(16) will be examined using simple physical arguments, which give results 
which are identical to those obtained in appendix A from the exact solution 
to the Poisson-Boltzmann equation of the second article in Ref. [20] . In order 
to accomplish this, let us write (p(z) as <j>(z) = 0o + <f>'{z) where 0o = 4>(D/2) 
and 4>'{z) is zero at z=D/2. Then, we can make Eq. (27) look like Eq. (16), 
if we identify uq with n s e~^°. Then we can write Eq. (27) as 



d 



% = -47r£ B [n (e-*' - (n s /n )e + *') - fcf>(z)). (28) 



dz 

When n s /n = « 1, Eq. (16) is definitely a good approximation to 
the problem, and we are justified in treating the system as one without 
excess salt. To determine the conditions under which Eq. (16) is a good 
approximation, we solve Eq. (16), in order to determine hq as described 
earlier in this section, and determine 0o from uq = n s e~^°, and use <po and 
the solution of Eq. (16) for <f>(z), which we identify with <f)'(z) (which is the 
approximate solution to Eq. (28)) to determine the conditions under which 
we may neglect e^ z > compared to e~^ z \ (Remember that <j){z) is negative.) 
Since it is easily seen from Eq. (10) that ^ is significantly greater than 
1 over much of the range of z from to h for high density polymer brushes, 
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all that is required in order to neglect is that n s /n be of order unity, 
which is already satisfied for the 0.1 M salt concentration (or 0.6 x lO 26 ??!^ 1 ) 
typical of living matter. It is demonstrated in appendix A that these results 
follow directly from the exact solution of the Poisson-Boltzmann equation 
with excess salt present [20J. 

6 Conclusions 

It has been shown using a modified version of the mean field theory of 
Miklavic Marcelja[19] for polyelectrolyte polymer brushes, which uses the 
non-linear Poisson-Boltzmann equation, that it should be possible for os- 
motic pressure due to the counterions to support a reasonably large load 
(about 10 6 Pa) with the tops of the brushes sufficiently far apart to pre- 
vent entanglement of polymers belonging to the two brushes, which has been 
argued to account for most of the friction. This load carrying ability is ar- 
gued to persist in the presence of an amount of added salt comparable to 
that found in living matter. Significant additional salt, however, provides 
screening which reduces the load carrying ability of polyelectrolyte brushes. 
Using counterions of higher valence will also not improve the load carrying 
because it will actually reduce the net charge on the polymers, by causing 
more counterions to condense [30J. This will in turn reduce the counterion 
concentration in solution. The load carrying ability of the brushes could be 
improved by using a solvent with a higher dielectric constant, which would 
reduce the value of the Bjerrum length, which appears in the denominator of 
the expression for the counterion osmotic pressure [Eq. (15)]. Using denser 
brushes or better solvents, which increase the ratio of brush height to polymer 
radius of gyration would improve the load carrying ability, since the mini- 
mum thickness of the interface region between the polymer brushes which 
avoids entanglements that lead to static friction, is given approximately by 
the quantity £ under Eq. (8) is proportional to h^ 1 ^ 3 . Making this region 
less thick increases the counterion osmotic pressure, as Eq. (15) shows that 
it is inversely proportional to the square of the interface region thickness. 
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7 Appendix A: The high and Low Salt Con- 
centration Limits of the Exact Solution of 
the Poisson-Boltzmann Equation 

The exact solution to the Poisson-Boltzmann equation [i.e., Eq. (20)] given 
in the second article in Ref. [20] for the electrostatic potential in the interface 
region between two polymer brushes in terms of elliptic functions, which are 
close to each other but not touching, is 

««>-^ wSw).») '- {1A) 

where z = D/2 — h and K = (87r£s^j) 1,/2 , where rij is the salt concentration, 
c n (x, k) is an elliptic function and k, the standard elliptic function parameter 
k [36], is given by k = [cosh(0(_D/2)/2)]~ 1 . First, let us consider the limit as 
<p(D/2) approaches zero. From Eq. (1A) we find that 

n{D/2) — (f)(D/2)cosh[K(D/2 — z)] (2A) 



cn(Kz, 1) 

since cn(x,k = 1) = sech(x). This is the solution in the Debye-Huckel 
approximation. 

Now let us consider the limit as \4>(D/2)\ approaches infinity. Since 
cosh(4>(D I '2) approaches infinity, k in Eq. (1A) approaches zero. It is easily 
shown that cn(u,k=0)=cos(u). Then, from Eq. (1A) we obtain 

sinh(4>(D/2)/2) 



sinh((f)/2) 



cn(Kzcosh((j)(D /2)/2),0) 
sinh((f>(D /2)/2) 



cos(Kzcosh{(j){D/2)/2)y 
which becomes for |0(£)/2)| >> 1, 

4>(z) w 0OD/2) + ln\cos 2 [(K(z - D/2)cosh(<j>(D/2)/2)]\. (3 A) 



Since as \<f)(D/2) \ approaches infinity, cosh(<ft(D /2)/2) is approximately equal 
to (l/_2)e^ D /V\/ 2 , Kcosh(<p(D /2)/2) becomes [27r£ B n ] 1/2 , with n = n s e^ D ^ 
n s e~^ D l 2 \ The latter result follows from the fact that since the counterions 
must have lower potential energy between the plates than outside the plates, 
4>(D/2) < 0. 
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